Exact free energy distribution function of a randomly forced directed polymer 
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We study the elastic (1 + l)-dimensional string subject to a random gaussian potential on 
scales smaller than the correlation radius of the disorder potential (Larkin problem). We present an 
exact calculation of the probability function V [F {u, L)] for the free energy F of a string starting at 
(0,0) and ending at (u,L). The function Vi^F) is strongly asymmetric, with the left tail decaying 
exponentially (ln'P(F — oo) cx F) and the right tail vanishing as \nV{F —>■ +cx)) oc —F'^. Our 
analysis defines a strategy for future attacks on this class of problems. 
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Over the recent years, the interest in static and dy- 
namic aspects of random manifold problems has in- 
creased significantly, with numerous applications in the 
field of random magnets vortex matter in type II 
superconductors dislocations in metals Q, charge- 
density waves in solids and many more. Also, the 
physics of random elastic manifolds is related to other 
fascinating topics, e.g., Burgers turbulence |6|, stochas- 
tic growth 1^] , or systems of interacting bosons [|| , all of 
which have attracted considerable interest recently. In 
this context, the Larkin model |^ has become a generic 
testing-field for various ideas, similar to the harmonic 
oscillator in quantum mechanics — both model systems 
are simple and exactly solvable, while at the same time 
being very instructive. Quite astonishingly, the study of 
the Larkin model has not been carried through in all its 
facets. In particular, the distribution fimction ViF) of 
the free energy F, an important element in the charac- 
terization of the dirty elastic manifold, is not known to 
date. The present work eliminates this deficiency with 
an exact calculation of P{F). 

A d-dimensional dirty elastic manifold (elasticity e) 
with n transverse degrees of freedom u(z) (u e M", z e 
R'') is described by the Hamiltonian 



H[u(z)] 



du 
dz 



C/(u(z) 



(1) 



The disorder potential C/ is a random gaussian variable 
with a correlator i^duj) decaying on the scale R^ , 



(L/(u,z)J7(u',z')) = K{\u-W\)5{z - z') 



(2) 



The partition function Z{u, L) of a string (d = 1) begin- 
ning and ending at the points (0, 0) and (u, L) can be 
written as a path integral 



(u,L) 



Z(u,L)= j V[n\z')] e^pi^-^J^dz' 



e I du! 



(3) 



(0,0) 



+C/(u'(z'),z') 



from which we obtain the free energy F = —Tin ^(u, L). 
Due to the random nature of the potential ?7, the free en- 
ergy _F is a stochastic variable and we are interested in 
its distribution function V(F), see Fig. 1. 

The imaginary time Schrodinger equation TdzZ — 
(T^/2e)dlZ - UZ, z e [0,L], satisfied by the partition 
function Z, directly maps to the Kardar-Parisi-Zhang 
equation Q d,F = {T/2e)AF - (VF)^ /2e+U describ- 
ing the stochastic growth of the 'surface' F(u, z), and 
differentiating the KPZ-equation with respect to u and 
going over to the gradient field v = VF leads to the 
Burgers equation [^. 
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FIG. 1. The function P(F) for the Larkin model (we 
have subtracted the trivial thermal component Fth from F; 
the parameters are /3 — e^u^/aL'^ and Fo = aL^/e, with 
a = —K"{0) the disorder parameter in the force correlator): 
the general case with /3 7^ ^ it 7^ (dotted line) is char- 
acterized by the tails P{F — + ~oo) ^ exp(— 7r^|_F|/Fo) and 
P{F ^ CX)) ~ exp{-l6F^/27P'^F§). For (3 = ^ u = 
(solid line) the distribution function vanishes identically for 
F > with P{F 0") ~ exp(-Fo/16|F|). The distribu- 
tion function for the random polymer problem on short scales 
u < Rk, see Eqs. (1) and (2), is independent of u and is given 
by the result for the Larkin model with (3 = 0. 



In studying random elastic manifolds, we are interested 
in the statistical behavior of the displacement field u and 
the free energy F as functions of the internal coordinates 



1 



z. Important results for the directed polymer problem 
concern the disorder induced line wandering (u^) cx L^^ 
and the free energy average {F) and its fluctuations. A 
more ambitious goal is the calculation of the probability 
distributions V{u, L) and T'{F, L). Only few rigorous re- 
sults are known, among them the behavior of (u^) and 
(F) of the Larkin problem ||ri|, the solution of the (1 + 1)- 
directed polymer problem by Huse, Henley, and Fisher 
I p^ , where C = 2/3 is an exact result, and the replica 
variational calculation for the [d + n)-random manifold 
by Mezard and Parisi ||l^ , who find ( = {4 — d)/{4 + n) 
for d > 2 and ( — {2 — d)/2 for c? < 2, an exact result in 
the limit n — > oo. The free energy probability distribu- 
tion function 'P{F) has been studied intensively (see Refs. 
|p^-|l8|) but few solid results are available: the left tail 
In-P(F -oo) cx -(|i^|VL)i/2 for the (1 + l)-problem 
(from which the free energy fluctuation exponent x = 1/3 
follows trivially) has been found by Kardar |l4| and some 
(partly contradicting) results for the tails of the velocity- 
difference correlation function of the Burgers equation 
have been reported recently |l^ . 

In this situation it is appropriate to look out for those 
problems which can teach us some definitive non-trivial 
results. Here, we investigate the exactly solvable model 
of a random directed polymer on scales smaller than the 
correlation radius of the random potential. The stud- 
ies on this model go back to the pioneering work of Larkin 
1^, who, exploiting the smallness of the displacement 
field u, expanded the potential in the Taylor series 



C/(u,z) w [/(0,z) - f(^) 



(4) 



This model is linear in u and thus allows for an exact 
analytical solution. Though simple, it produces a num- 
ber of non-trivial results, of which the distribution func- 
tion P{F) is still unknown. Below, we solve the problem 
through a mapping on a system of interacting bosons. 

We start from the partition function (^ for the 
(1 + l)-dimensional string in a random force field f{z) 
|P0| , with f{z) a random gaussian variable with cor- 
relator {f{z)f{z')) = a5[z - z'), a = -K"{0), see 
Eqs. (|[) and (^) [the generahzation to the (1 + n)- 
problem is trivial]. Replicating the partition function n 
times wc obtain the Green function G{ui, . . . ,Un]L) = 
{Z{ui, L) . . . Z{un, L)) satisfying the imaginary time 
Schrodinger equation (z e [0, L]) 



dz 2e ^ duf 

-1 



(5) 



with the initial condition G{ui, . . . , Un] z ^ +0) — 5{ui) 
S{u2) ■ ■ ■ 6{un). The boundary conditions at infinity read 
G(ui, . . . , — > oo,...,M„;z) 0. Equation (j|) has 
the form TdzG = —HG, with H the Hamiltonian de- 
scribing the imaginary time evolution of a system of in- 
teracting bosons, with the temperature T playing the 



role of Planck's constant and the internal coordinate 
z of the string mapping on to the imaginary time co- 
ordinate of the boson system. In the new variables 
t = ^al2tTz and Xi = {ae/2T^Y/'^ u„ Eq. (||) takes the 
simple form dtG = (1/2) Y,i dfG + J2i,j ^i^j G. The spe- 
cific form of the Hamiltonian motivates an Ansatz G — 
g{t) exp[(l/2) aij{t)xiXj] for the Green function, re- 
ducing the problem to a system of ordinary differential 
equations g/g = (1/2) Yjk '^kk and dy = o-ikajk + 2. 
We assume the matrix A = {oij} to be replica symmetric, 
i.e., an = a and Ui^j = b (this assumption implies no re- 
striction, as the present problem has a unique solution). 
The three resulting differential equations g/g = na/2, 
b = {n- 2)6^ + 2ab + 2, and d = {n - 1)6^ -1-0^-1-2 are 
easily solved: concentrating first on the difference a — 6 
we find, using the initial condition at i = 0, the result 
a — b — —1/t. Using the Ansatz c{t) — exp[— n J dt' b{t')] 
in the equation for 6, we can reduce the problem to 
the Bessel equation and find the solution b{t) = 1/nt — 
2tcot{^/2nt)/^/2rit, where we have again made use of 
the initial condition. Finally, the result for g(t) reads 

g{t) = (2n)i/4yiexp[-(l/4)lnsin2(V2^0] and assem- 
bling the various elements, we arrive at the following final 
expression for the moments Vin) — (Z"(2;, t)/Z^(a;, i)) 
of the replicated partition function (the arguments in V 
are x u and t L; we subtract the free energy Fth of 
the thermal model (U = 0) as we are only interested in 
the difference AF = F — Fth with and without disorder; 
in the following we drop the symbol 'A'), 



V{n) = exp 



Tix^ 1 , sm (\/2nt) 

In 7^ 

2t 4 2nt^ 



9 9 cot(v2nt) 



2nt 



(6) 

The moments (^) give us access to all characteris- 
tic quantities of the random directed polymer on short 
scales. In particular, the n — > limit of (^) determines 
the free energy average (F) [|n],| and the probability 
function V{u,L) for the displacement field u The 
large n limit determines the 'tails' {F ±oo) of the dis- 
tribution function V{F)] the 'left tails' V{F —oo) for 
the random directed polymer problem have been stud- 
ied on short 0,11^ and long scales |]l^-|l6|. However, in 
order to determine the 'body' of the function 7'(F) one 
has to calculate all the moments (||) (in the boson lan- 
guage the determination of the 'tails' and the 'body' of 
the distribution function V{F) amount to calculating the 
ground state of the n-boson problem and the full spec- 
trum, respectively). In the following, we present a brief 
derivation of these quantities as they follow from 

The n —> limit of (|^) determines the disorder aver- 
aged free energy [|| (F) = -T(\n[Z/Zo]) = -Tlim„^o 
[(Z/Zo)" - l]/n = -aLVl2e. Similarly, this limit deter- 
mines the correlation functions for the displacement field 
u: These follow from the probability function V{u, L) = 



{Z{u,L)/j:^duZ{u,L)) 



-^v?), where 



2 



7~"^ = 2(TL/e + aI/'^/3e^). The mean squared displace- 
ment (tt^(I/)) = {[u{L) — u(0)]^) on scale L then is given 
by {v?{L)) = {T/e)L + {a/'ie^)L'^ and crosses over from 
diffusive thermal- to random wandering at a length scale 
L ^ (T I a. 

We turn to the calculation of the probability distri- 
bution function V{F). The moments V{n), see (^), are 
nothing but the Mellin transform of V{F) pft], 



V{n) 



+00 



dFV{F)eM-nF/T), (7) 



and we have reduced the problem of calculating V{F) to 
that of inverting the Mellin transform. The latter de- 
pends only on the combinations A = 2t^ n = (aL'^/eT) n 
and P = /2t^ — e^v? /aL^ ^ i.e., we can rewrite (|^) in 
the form [Fo = 2Tt^ = aL'^/e) 

/+00 
dFV{F)eM-^F/Fo) {^V{n)) (8) 
-00 



exp 



1, sm'iVX) B^col{^/\) 

—A — 7 m —A ^ — 

2 4 A 2 VA 



The function /(/3, A) is analytic in the halfplane Re{A} < 
TT^ with /(/3, A — > TT^ — 0) 00. We then can invert the 
Mellin transform and write the function ^(F) as an in- 
tegral of / over the imaginary axis, once we have found 
the proper analytical continuation of A) (see below). 



V{F) = 



1 



2^Fo 



dA/(/3,^A)exp(^AF/Fo). (9) 



Note that the condition for the existence of the inverse 
Mellin transform, J d\\f{(3,iX)\ < co \s satisfied: for 
/3 = we have |/(0,iA ±i(x,)\ - exp(-^|A|/2) while 
for /3 ^ 0, |/(/3,iA ^ ±i^)\ ~ exp(-/?|A/2|3/2). 

Tails: With A) analytic for Rc{A} < tt^ the ex- 
ponential exp(— AF/Fo), F < 0, in (||) is compensated 
by the decaying tail of 'P(F) up to the point A = tt^, 
implying a left tail of the form 



V{F -00) - exp(-7r2|F|/Fo). 



(10) 



The right tail of 7'(F) is determined by /(/3, A —00). 
Since /(/3 7^ 0, A -> -00) ~ exp(/3|A|3/2/2) remains fi- 
nite at finite |A|, we conclude that V{F) decays faster 
than exponential. We can invert the Mellin transform 
asymptotically via the method of steepest descents and 
obtain 

V{F -> c») - exp(-16FV27^2^o^) ^ 0). (11) 

The case /3 = has to be treated separately: With f{(3 — 
0,A ^ -cx)) - exp(-vlA|/2), we find V{F > 0) ee 0. 
Using the Ansatz P{F) oc exp(— A/|F|") in combination 
with the method of steepest descents produces the result 



p(F^-0)-cxp(-Fo/16|F|) (/3 = 0). 



(12) 



Body: We have to find the proper analytic continua- 
tion of the Mellin transform ^ to imaginary A values 
(imaginary boson number), see (^). Inspiration how to 
carry out this analytic continuation can be obtained via 
the alternative route defining the distribution function in 
terms of a path integral over the stochastic field /, 

V{F) - J V[f{z)] S{F - n[uo{.f)] + e,?/2L), (13) 

where uq{z) is the solution of the equation euQ(z) ~ 
— f{z) with the appropriate boundary conditions Mo(0) = 
and uo(F) — u and 7i[ito(z)] is the associated energy 
{{x,y) denotes the usual scalar product). 



-n[u,{z)] = 1^ + ^(/,i/) - l{z, /), (14) 



Af 



dz' 



L- z 
L 



(L — z') -\- max(z, z') — L 



We determine the spectrum of the Hermitian operator A 
{An = —L'^/ir'^n^) and after integration over the stochas- 
tic field /, we arrive at the distribution function in the 
form V{F) = /+^(d/i/2^) exp[i^F - g{p)l with 



In 



aLu II 



2, ,2 r 



n -1' 



(defining = iX/Fo, a few algebraic manipulations trans- 
form this result back to (||)). The above result allows 
for a straightforward inversion of the Mellin transform 
via numerical integration and thus to reconstruct the full 
distribution function 'P{F), see Fig. 1. Furthermore, it 
implies that 7'(F > 0) = for u = 0: We integrate over 
the contour in the //-plane made from the real axis and 
the upper (F > 0) semi-circle. The integral over the arc 
vanishes asymptotically. With all branching points sit- 
uated in the lower half-plane, the contour integral van- 
ishes and 7'(F > 0) = 0. This property of P(F) can 
be understood in terms of a minimization of the func- 
tional 'H[u{z)]: With the endpoint (0, L), the test func- 
tion ut{z) = satisfies the boundary conditions and we 
have H[uo] < H[ut] = 0, hence 7'(F > 0) = 0. This anal- 
ysis cannot be applied to other endpoints with u =/= 0. 

The distribution function V{F) for the Larkin model, 
particularly its 'tails', have been studied before and a 
comparative discussion is instructive. Expanding the 
correlator K{u) in a Taylor series, {U{u, z)U{u' , z')) « 
[^(0) - (a/2)(u - u'Y]5[z - z'), the authors of Rcfs. 



|17| , |18[ arrive at the asymptotic expression (Z") ^ 
exp[— C\/n(n — 1)F], which cannot be identified with 
the Mellin transform of some distribution function 7'(F), 
however. Keeping only the second term exp(C\/ni), as 
suggested in Ref. |l8j, is questionable: though the left 
tail ~ exp(— i?L2/|F|) reproduces the desired answers 



3 



for the fluctuations of the free energy SF oc and the 
hne wandering u cx L^/^, this tail is unacceptable for a 
properly normalizable distribution function. The origin 
of the problem seems to lie in the correlator expansion 
itself: determining the exact Mellin transform using the 
expanded correlator we find the result 



cxp 



(n) — exp 



1 



In 



K 



4 sinh^ {y/2^t) 2 



where K = K{0)L/T^. While the leading term has the 
correct sign, the fact that Pexpin = — 7r^/2t) = im- 
plies that VexpiF) assumes negative values, in conflict 
with the required positivity of V. The findings of Parisi 
pTf show that the expansion provides correct results for 
the mean value of the free energy and the displacement 
correlation function, however, higher moments emphasiz- 
ing larger distances u — u' , where the expanded correla- 
tor becomes negative, turn out wrong and so does the 
function V{F). Carrying the expansion of the potential 
further, U{u, z) « t7(0, z) - f{z)u + g{z)u^ /2, the cal- 
culation of the distribution function can be done 
along the lines described above (see (|l^)) and we recover 
the result for the Larkin model for /3 = 0, see Fig. 1: as 
expected, the answer for the translation invariant ran- 
dom polymer problem and (|^) with a homogeneous 
correlator K{\u — u'\) is independent of the coordinate u 
of the polymer's endpoint. 

Next, we address the large scale behavior of the ran- 
dom directed polymer, which can be solved via the Bethe 
Ansatz technique . The ground state energy Eq (n) cx 
produces the left tail lnV{F) cx -|F|3/2/ii/2 
and provides the scaling SF cx L^^"^ for the fluctuations 
in F and the correct exponent C = 2/3. Still, the gen- 
eralization of this procedure seems not straightforward. 
E.g., for the (1 -I- 2)-problem the ground state energy 
increases faster than any power of n, implying that SF 
behaves logarithmically and ( — 1/2, in contradiction 
with the expected superdiffusivc (C > 1/2) behavior of 
the (1 -|- 2)-dimensional string. It then seems important 
to know the complete Mellin transform 7'(n) for an ac- 
curate calculation of the wandering exponent 

An interesting remark concerns the wandering expo- 
nent ( ~ 3/2 for the Larkin problem and its relation 
with the Kolmogorov scaling for the small-distance ve- 
locity fluctuations in the Burgers turbulence problem, 
{[v{u,L — oo) — v{—u,L ~ oo)]^) cx u^'^, with u = 1/3 
(see also [^). The consistency between the two exponents 
( and is easily checked: The fields in the two prob- 
lems are related via v = 9„F ~ F/u, hence F cx u^^'^. 
In the Larkin model F (x (x w^/^ and we obtain 
2/C = 1 + ly, which is indeed satisfied. However, this 
result does not carry a deep physical meaning but is sim- 
ply a consequence of scaling. Physically, the two prob- 
lems are very different: In the Larkin model we study the 
free energy F both as a function of u and (finite) L with 



only one pinning valley relevant, whereas for the case of 
Burgers turbulence the stationary correlation functions 
mk{u) = {[v{u,L = oo) — v{—u,L = oo)]*') are studied, 
where many disorder-induced valleys are relevant (in the 
language of the directed polymer problem). 

In the end, the calculation of the probability distribu- 
tion V[F{u, L)] turns out to be quite a non-trivial prob- 
lem: a successful attack on the problem requires the pre- 
cise knowledge of all the moments V(ri). Under these 
circumstances, we then can proceed with the analytical 
continuation to imaginary boson number and invert the 
Mellin transform to obtain P{F). The successful appli- 
cation of this strategy to the Larkin model uncovers the 
limitations of previous attacks on the problem and may 
serve as a guideline for future studies. 

We thank Thomas Nattermann and Valerii Vinokur for 
helpful discussions. 
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We omit the term U{0, z) dz corresponding to the triv- 
ial part of disorder. The Mellin transform of the proba- 
bility distribution for the full free energy is given by the 
product V{n)exp{K{0)Ln^ /2T'^). 

The Mellin transform /(p) of the function f{x) is de- 
fined as /(p) — j/^ dx f{x)x^~^ . With x — exp(— f) this 
transforms to /(p) = Jj*""" /(exp(— i)) exp(— pt). 
Cutting off the growth of u at Rk, the result crosses over 
to the extensive behavior {F{L > Lc)) ~ Uc{L / Lc) , with 
Lc = {t^R\/af'^ and Uc = {eaR%Y^^ (^ collective 
pinning theory, see Ref. jpl). 
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